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TECHNICAL NOTE NO. 1809 


THE METHOD OF CHARACTERISTICS FOR THE DETERMINATION OF 
SUPERSONIC FLOW OVER BODIES OF REVOLUTION 
AT SMALL ANGLES OF ATTACK 
By Antonio Ferri 

SUMMARY 


The method of characteristics has been applied for the determination 
of the supersonic— flow properties around bodies of revolution at a small 
angle of attack. The system developed considers the effect of the vari- 
ation of entropy due to the curved shock and determines a flow that 
exactly satisfies the boundary conditions in the limits of the simpli- 
fications assumed. Two practical methods for numerical calculations are 
given. 


INTRODUCTION 


For the determination of aerodynamic properties of bodies of revo- 
lution at supersonic speeds, two methods have been used: a method that 

uses the small— disturbances theory and a method that uses the charac- 
teristics theory. Both methods are successful in the determination of 
the flow properties for bodies at zero angle of attack, but the precision 
of the small-disturbances theory decreases when a body of revolution at 
an angle of attack is considered. 


For bodies of revolution having supersonic flow everywhere, the 
theory of characteristics can also be used at an angle of attack. 


The method of characteristics for the determination of the flow 
field around bodies of revolution at an angle of attack was first used 
by Ferrari (reference 1) in 1936. Ferrari considers the flow as potential 
flow and develops a method for the analysis of the flow field around a 
body that in the approximation of potential flow appears to be general 
and can be applied to bodies of any shape and with any angle of attack. 

In the determination of the flow properties along the first character- 
istic surface from which the analysis starts, however, Ferrari analyzes 
the flow around a cone of revolution; and in this part of the analysis 
only small values of angle of attack are considered. 

Sauer in 1942 (reference 2) considers the same problems and shows 
that for small values of angle of attack, the analysis of the flow 
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field around a body of revolution can be made by applying the charac- 
teristics method only in one meridian plane; and, therefore, Sauer uses 
characteristic lines in place of the characteristic surfaces considered 
by Ferrari. Sauer, in the development of his system, is interested 
essentially in the analysis of the flow around circular cones; and when 
the method is applied to bodies of revolution of shapes different from 
cones, the boundary conditions are no longer satisfied. The flow obtained 
from the solution used, also at small angles of attack, wets a body that 
is not a body of revolution. The body can be obtained from the body of 
revolution considered initially by curving its axis of symmetry. Sauer 
also assumes that the flow is potential flow. With this assumption, the 
flow must be considered as potential flow for the case of the body at 
zero angle of attack also; therefore, all the effects of entropy gradients 
are neglected. 

The flow field around circular cones at small angles of attack has 
been analyzed in a more exact form by Stone. (See references 3 and k.) 

In his analysis. Stone considers the flow as rotational flow and, there- 
fore, takes into account the effect of entropy gradients on the velocity 
distribution. This effect exists only when the cone has an angle of 
attack and, at low Mach numbers, is small but of the same order as the 
effect of other parameters that are considered in the analysis. 

Here, the method of characteristics is extended to the analysis of 
the flow field around a body of revolution at small angles of attack for 
the case of rotational flow. The effect of entropy gradients about bodies 
of revolution even at small angles of attack can be important because the 
entropy gradients that exist in the stream for small angles of attack are 
due to the variation of curvature of the shock existing at zero angle of 
attack also, together with the fact that the shock surface does not have 
axial symmetry with respect to the direction of the undisturbed velocity. 

The method presented permits the determination of a flow that in 
the assumption of small angles of attack exactly satisfies the boundary 
conditions and, therefore, wets the body of revolution considered. This 
method is given in a form that permits its application to practical 
problems and requires either numerical or numerical and graphical calcu- 
lations of the same type as the calculations used for the analysis of the 
flow around bodies at zero angle of attack. The method can be applied 
to cases in which the entropy variations can be neglected or are zero. 

In these cases the terms that contain the entropy variations become zero. 


SYMBOLS 


x,y,9 cylindrical coordinates (fig. 1) 

r,\|f,9 polar coordinates (fig. 8) 

7 local velocity (function of x, y, 9) 
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u,v,w 


v r ^ n ,w 


P 

P 

7 

a 

a 

3 

<P 

T| 


a 




N 


Vip 


w 


a 

AS 

n 

N 


velocity components in cylindrical coordinates (u along 
x-axis, v along y-axis, and w normal to meridian 
plane ) 

velocity components in polar coordinates (v r along r, 
v n normal to r in meridian plane, and w normal to 
meridian plane) 

limiting velocity corresponding to adiabatic expansion to 
zero pressure 

pressure 

density 

ratio of specific heats 

speed of sound ^a 2 = 7p^ 

angle of attack of body 

Mach angle ^sin p = ^ 

angle between velocity V and x-axis 

angle between the axis of the cone tangent to the shock 
and the axis of the body 

angle at the apex of the cone tangent to the shock 

tangents to the characteristic surfaces in the meridian 
plane 6 = Constant 

velocity component normal to the shock surface 


velocity component along the generatrix of the cone 
tangent to the shock 

velocity component tangent to the cross section of the 
cone tangent to the shock 

angle between the tangent to the shock and the axis of 
the body 

entropy variation for unit mass 

normal to the streamline in the plane 6 = Constant 

normal to the surface of the shock 
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H,L,K,Z 

A 

Q 

Di^D2 


coefficients defined "by equations (24) 
coefficient defined by equation (4l) 
coefficient defined by equation (45) 
coefficients defined by equations (55) 


Aj_, A2 ^Pi^P 2^T coefficients defined by equations (60) 


E radius of the holograph diagram 

Subscripts: 

0 free— stream flow quantities 

1 flow quantities for the condition of zero angle of 

attack 

2 flow quantities related to the effect of angle of 

attack as defined .in equations (5), (6), and (8) 

The prime ( 1 ) represents quantities in front of the shock and the 
double prime ( ft ) represents quantities behind the shock. 


EQUATION OF MOTION FOE FLOW AROUND A BODY OF 
EE70LUTI0N AT A SMALL ANGLE OF ATTACK 


Consider a cylindrical coordinate system in which the x— axis is 
coincident with the axis of the body of revolution, the y— axis is normal 
to the x--axis in any meridian plane, and the position of every meridian 
plane is defined by the angle 6 measured with respect to the meridian 
plane that contains the direction of the undisturbed velocity (fig. 1). 


Euler f s equations of motion for steady flow in cylindrical coordi- 
nates are: 


1 dp _ bu. 
p bx bx 


u + 


8u 


v + 


8u 

r be 
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(la) 


1 by _ 8v 
p by 8x 
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8v 

dy 


V + 


8v 

r bG 
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w^ 

y 


(ib) 


1 by 
P 7 bO 
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8x 


u + 


bw 

8y 


V + 



w + 


vw 
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The continuity equation in cylindrical coordinates can be expressed 
in the form 


d(pu) + S(pvy) + d(wp) = Q 
bx y dy y b9 


whereas the law of conservation 


of energy 


can be written in the form 




<*£ _ _E_ 

Bx p 2 Bx/ 




+ w 



(3a) 




& _ JL 
3y p 2 byj 


u 


+ v ^ 
ay by 


+ w 



(3*>) 


7 /I dp p_ dp \ 

7-1 \p y be p 2 y bej 



+ V 


dv dw \ 

+ w I 

y bd y bd/ 


(3c) 


If the density is eliminated from equation (2) by means of 
equations (1) and (3) and the quantity a is introduced defined by 



the following equation can be obtained: 


du 

dx 




V 

y 


vu / du , dv\ 
a? \by bxj 


uw /dw du \ _ vw / dv dw\ _ Q 
a? \dx + y bej a 2 \y d0 byj 


00 


In this analysis only small angles of attack will be considered, 
and, therefore, only the first— order effect of the angle of attack will 
be determined; whereas the quantities of the same or higher order than 
the square of the angle of attack will be neglected. In this approxi- 
mation the velocity components of the flow around the body can be 
expressed in the following form (references 1 to k) : 


u = uq + au2 cos 9 (5a) 

v = v^ + 1 x ^2 cos 9 (5*) 
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w = 0CW2 sin 9 


( 5 c) 


where u, v, and w are functions of the three coordinates 1, y, and 9; 
whereas u^, U2, ^ 2 .’ an< ^ w 2 functions only of the coordinates x 

and y of any meridian plane. The quantity a is the angle of attack 
of the body, the quantities with subscript 1 are the quantities 
existing at the position (x,y) for the body considered at zero angle of 
attack at the same Mach number, and the quantities with subscript 2 are 
functions that take into account the effect of angle of attack. 

It will be shown in the following considerations that the form 
assumed in equations (5) for the velocity components permits the boundary 
conditions to be satisfied in the simplifications assumed. From 
equations (5), for small angles of attack. 


whereas equation ( 4 ) becomes 



Equation (7) is similar in form to the corresponding equation for 
the case of the body at zero angle of attack and differs only in the 
term — . In order to analyze the differences between this expression 

and the expression for the axial symmetrical case and in order to obtain 
another relation that defines the quantity w, the relation between 
rotation of the flow and entropy gradient will now be introduced. 

For perfect flow the entropy variation at any point can be expressed 
in the form 


For small angles of attack, therefore, by the use of equations (6), 


p = Pi + ap2 cos 9 


(6a) 


p = + ap2 cos 9 


(6b) 
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1 -1- a 52 cos 9 — 7a — cos 9 


Pi Pi 
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or 


and 


AS = AS^ + a AS 2 cos 0 


ASg = ( El - 

r-1 V Pi Pi y 


( 8 a) 


(8b) 


where AS^ and ASg are functions only of x and y. 

Between rotation of flow and entropy the following relation exists: 


curl V x V = — grad S 
7 R 


or for small angles of attack 


8 S a^ _ _ y / 8 v _ chA 
8 x 7 R 8 y/ 


(9a) 


8S a£ _ / 8 v _ 8 u\ 

8 y 7 R \ 8 x 8 y/ 


(9b) 


c)S a 2 _ v / 8 w _ _ 8 v_\ + wv _ / 8 u _ 8 w\ 
y 89 7 R \ 8 y y 8 9/ y \y 89 8 x/ 


(9c) 


If n is the normal, in the meridian plane 9 = Constant, to the local 
tangent to the streamline, then 


8 S _ 8 S t t 8 S u 
8 n 8 x V 8 y V 

= 2 E y _ 8 u\ 

a 2 \ 8 x 8 yy 


( 10 ) 
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whereas from equation (9c), when equations (5) and (8) are used, it 
follows that 


=£ = + I2\ + + /U£ + bwr\ 

7K 7 \dy jJ y \y dx/ 


(ID 


or 


_ 8w2 , „ c>W2 _ VW2 + VT2 + UU2 a? A „ 

v — — + u -r — — — — Abo 

oy dx y y 7 E * 

Equation (7) can he written in the following form: 


( 12 ) 


du A _ uM , d 7 /1 v 2 ) , v _ 2ut dv . uv dS , dw _ 

3xV e?) a 7 V s?) 7 a*&c (13 

Equation (13) together with equation (12) defines the law of motion of 
the flow around the body at small angles of attack. This equation will 
he used as a basis for the calculation of the flow field by the method 
of characteristics to be treated in a later section. 


CONDITIONS AT THE SHOCK FRONT 


Equations (5) and (6) represent a stream that wets a body of revo- 
lution at a small angle of attack. In order to satisfy the boundary 
conditions at the surface of the body, the functions U2* V2, and W2 
must be properly selected. Equations (5) and (6) must, however, satisfy 
the boundary conditons at the shock surface also in order to be a 
solution of the problem. It is necessary, therefore, to show that a 
shock surface can exist across which the undisturbed stream inclined 
at a with respect to the axis of the body is transformed into a flow 
represented by equations (5) and (6). 

In order to show that the shock boundary conditions can be satis- 
fied, the following procedure will be employed. A shock surface distorted 
in a manner to be described is assumed. Then, the free— stream velocity 
ahead of the shock will be resolved into three components: v^ normal to 

the shock, Vrp tangent to the shock in the plane 9 = Constant, and w 
perpendicular to the plane 9 = Constant. Similarly, the flow behind 
the shock will be resolved into three components. In addition, each 
component of the flow behind the shock will be divided into two terms: 
one term for zero angle of attack and one term for the difference due 
to the angle of attack (for example, u = uj_ + au2 cos &)• Then, 
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the conditions of equilibrium at the shock will he imposed, and it will 
he shown that the terms u 2 , v 2 , and w 2 at the shock are independent 
of 9 when the angle of attack is small as initially assumed; hence, 
the distorted shock is consistent with the field of flow behind it. Such 
a shock surface can be obtained by deforming the shock surface produced 
by the body when the angle of attack a is zero in the following way 
(fig. 2): 

When a = 0, the shock surface is a surface of revolution in axis 
with the body; therefore, if OP, OP 1 is the curve intersection of the 
shock with the meridian plane 0=0, then for a = 0, the tangent AQ 
at any point Q of the curve OP is the generatrix of a circular cone 
having the vertex at a point A of the axis and tangent along the 
circle QQ ? to the shock surface. The shock surface, therefore, can be 
considered as a surface envelope of circular cones having the axis 
coincident with the axis of the body but having variable cone angle and 
variable position of the apex A along the axis AB of the body. For 
the case of a ^ 0 the shock surface is not a surface of revolution but 
can still be considered, for small angles of attack, as the envelope of 
the same circular cones considered for the case a = 0. These cones 
have the same apexes and the same cone angles as the cases for a = 0 
but do not have the axis of symmetry AB coincident with the axis of 
the body AB although they are rotated in the plane 9=0 with respect 
to the body axis. The angle tj through which each axis of the cones 
must rotate in the plane 0=0 with respect to the axis of the body, 
is not constant but varies for each cone considered. For example, the 
cone AQQ ? tangent to the shock surface for a = 0, when a ^ 0 3 must 
be rotated by an angle T) to the position AQpQp 1 ; the axis ABp 
remains in the plane 0=0. 

The shock surface so generated is consistent with the flow repre- 
sented by equations (5) and (6), and this can be shown in the following 
way: 


Consider a point P of the shock produced by the body at an angle 
of attack, and consider the cone tangent to the shock at the point P 
(fig. 3). Call a the semiangle of the cone with respect to its axis 
of symmetry. The axis of this cone is inclined at an angle q with 
respect to the axis of the body and lies in the plane 0=0. 

The uniform velocity V Q ahead of the shock is decomposed in the 
three components: v^ f in the direction PB normal to the shock, Vp 1 

in the direction AP along the generatrix of the cone, and w* in the 
direction normal to the plane APB. These components are, at small angles 
of attack. 


(l4a) 


V N* = V 0 sin cr — V 0 (a — t})cos a cos 9 
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v T * = V Q cos a + V 0 (a — n)sin a cos 9 (14b) 

W = -V 0 (a - r))sin 9 (l4c) 

Strictly, in equations (14) i|r must be written in place of 9; 
but, for small angles of attack in equations (14), the difference 
between \|r and 9 can be neglected. Indeed, 


t = 5 + 0* 


where 5 is of the order of a, and 0 » differs from 9 by a quantity 
of the order of a. 

The Telocity components behind the shock are (fig. 3) 

V = u sin a — v cos a + T| cos 9(v sin a + u cos cr) ( 15 &) 

Vrp" = u cos a + t sin a + T) cos 9(v cos a — u sin cr) ( 15 b) 

w" = w + T)(u — t cot cr) sin 9 (15c) 

where u, v, and w are the velocity components behind the shock in 

cylindrical coordinates in axis with the body at the point P considered. 

The velocity components u, v, and w at the point P can be 
expressed in the form given by equations ( 5 ), in which the components uq 
and Vq are the quantities obtained at the point P for the condition 

of a = 0 and are direct functions only of x and y. The point P, 

however, is a point of the shock, and its coordinates x and y change 
when the coordinate 9 changes; therefore, the velocity components uq 
and Vq at P also change with 0 q. In order to separate the part of 
the components u, v, and w dependent on 9 from the part independent 
of 9, the velocity components uq and vq at P will now be expressed 
as a function of the flow properties at a point Pq near P, having a 
constant value of x and y for every value of 9. 

Now, it has been assumed that the angle of the cone a tangent to 

the shock at the point P is equal to the angle of the cone tangent to 

the shock for the condition of zero angle of attack at the point Pq 
(fig. 4). The point Pq is obtained on the shock by rotating the 
cone APQ tangent to the shock for the condition a = 0 through an 
angle T| around the axis AN normal to the plane 9 = 0 at the apex A 
of the cone. Because for the condition of zero angle of attack the 
velocity components uq and vq are independent of the coordinate 9, 
the velocity components uq and vq at P(xp,yp) (fig. 4) are equal 
to the velocity components at I*2( x p^yp)j in the plane AP^C. Therefore, 

if AN is the distance P 2 Pq, 


u lp 




AN 


( 16 a) 
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where (fig. 4) 



X P 1 T 1 

AN = cos 0 

cos 0 


Substituting equations (5), (l6), and (17) in equations ( 15 ) 
in 


v N " = (uj_ sin 0 — vq cos 0)p^ 

+ a cos 0(u2 sin 0 — vp cos c)p^ 

+ q cos e(uq cos 0 + vq sin s)p 1 

■^Pl 1 ! /chiq ovq \ 

cos 0 ^N dN J r 1 

Vp" = (up cos 0 + vq sin 0 )p^ 

+ a cos e(u£ cos 0 + v 2 sin 0) p 

— t) cos e(up sin 0 — vp cos c)p 


fal l dv-, \ 

+ — - — COS 0 — — cos 0 + — 1 sin 0 p, 

\8n sn r * 


COS 0 


w" = aw 2 p sin 0 + q (uq — vq cot c)p sin 0 


( 1 6 b) 


(17) 

results 


(l8a) 


(l8b) 


(l8c) 
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For the condition of zero angle of attack at Pi_ 



= sin a — cos a 


( 19 a) 



cos a + sin a 


(I9h) 


and for the condition of the equilibrium at the shock at zero angle of 
attack 


Y n cos a = Tm " 


(20a) 


V Vo sin 0 = (vf _ 


(20b) 


At the point P for the case of a small angle of attack. 


(VV)p - ( 7 ! 2 - V 2 )j 


(21a) 



(21b) 


w 


T! 


W f 


(21c) 


MCA TN No. 1809 


13 


If equations (19) are used, equations (l8a) and (l8b) can be 
written in the form 


(V)p = ( v N 1 ")p 1 + a 003 9 ( v N 2 , ')p 1 + C0S 9 ( t Ti" 


pi 


X P x Tl 

COS 0 


cos 9\ 


% 


8 n /p 


( t t")p - (V)p + a 008 9 (V)p, - 1 006 9 (V)ri 


Xp T ) / OV T 

+ COS 9' ± 


cos a 


8N /p. 


Therefore, from equations (14), (20), and (21), 

/ xp q dv N "\ 

-(a - T))V 0 008 o(v Nl ") pi * T 0 sin <r(ar N2 " + n v Tl ” + 


_ _ £(r -Ll y : 2 C0S 0 B 0 ^ a — ^ ) 

7 + 1 


(22a) 


(a - q)V 0 sin a = (V) a - ^ ) Pi ( 22b ) 


or 


% 


o /Pi 


, t) \ . t] /y 

1 ] sin a + — 1 — i 


a 


1 x *i 1 ( 8 V' 


a \ V Q /■p 1 a sin 0 V Q \ SN / p^ 


(23a) 
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-A- -2(L^)(i 

tan a \7 + 1/ \ 



sin a + 





a 



a cos a V 0 



Pi 


(23b) 





— U]_ + cot 

v: 



a 

a 


(23c) 


In equations (23) the coordinate 9 does not appear; therefore, 
for the shock considered, the functions U2, V2, and W2 are independent 
of 9, and equations (5) represent a flow condition in agreement with the 
conditions at the shock. 

The ratio j]/a which appears in equations (23) is independent 
of a; therefore, for a given point P-j_, T]/a remains constant in all 
the range of angle of attack in which the simplifications assumed are valid, 

^Indeed, U2* and W2 are also independent of the angle of attack 
(equations (6)^. The values of u^ v^, w^, and t\/cl must therefore he 
determined only for one value of the angle of attack. 


METHOD OF CHARACTERISTICS FOR FLOW AROUND A BODY OF 
REVOLUTION AT A SMALL ANGLE OF ATTACK 


In this section the method of characteristics is applied to 
equation (13) to establish equations which will permit the flow field 
behind the shock to be calculated by a point— by— point process. If the 
flow is anywhere supersonic, equations (12) and (13) permit the determi- 
nation of the flow around a body of revolution at a small angle of 
attack by using the method of characteristics. Equation (13) can be 
written in the following forms 

h |u + L^v: + 2K^+Z = 0 

ox ay OX 
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Z _ uv _1_ dS + dw + v 
T 7E Sn y B0 y 


(24) 


If cp is the angle between the velocity V and the x— axis 
and p is the Mach angle , 


tan cp 


>/v^ + 

U 


sin P - y 


or for small angles of attack 


tan qp = — 
u 


and 


sin p = 


s/v 2 + u 2 
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The tangent to the line intersection of a characteristic surface 
with the meridian plane 9 = Constant is 


X a = tan(cp + 0) = - - i >/k 2 - HL 
H H 

X“b = tan(<p — p) = - + i s/k 2 - HL 
H H 

where X. a is the tangent to a line corresponding to the characteristic 
surface of the first family and Xq, is the tangent to a line corres- 
ponding to the characteristic surface of the second family. The 
terms X. a and Xq, are solutions of the equations (reference 5) 

HX 2 -2K\ + L = 0 (26) 

Because u, v, V , and a can he considered to he given by an equation 
of the type of equations (5)* 9 and 0 can also he written in the 
form 


(25) 


qp = cp]_ + ooqp 2 cos 9 
0 = 0^_ + a02 cos 9 


The characteristic surfaces are not, therefore, surfaces of revolution 
hut can he obtained, as was true for the case of the shock, as an enve- 
lope of circular cones with their apexes at the axis of the body and 
their axis of symmetry in the plane 9=0 and inclined with the axis of 
the body. 

The determination of the u and v components of the velocity in 
any point of the flow can he obtained hy equation (13) hy performing a 
transformation in order to obtain a law of variation along the character- 
istic lines (reference 5). Indeed, for every point of any meridian 
plane (for example, of the meridian plane 9=0, or 9 = ^) two charac- 
teristic lines can he obtained as the intersection of two characteristic 
surfaces with the meridian plane. Along these lines the variation of 
the u and v velocity components is determined hy the equations of 
characteristics that can he derived from equation (13)* Assume that in 
two points P^ and P 2 (fig* 5) of the meridian plane 9 = Constant 
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^for example, 9 = 0, or 9 = the velocity components are known. 

From equations (25) the tangents to the characteristic surfaces in this 
meridian plane can be drawn and the velocity components u and v at 
the point intersection of the two tangents, can be obtained in the 

first approximation. 

The equations of characteristics can be obtained by analyzing 
equation (13) along the characteristic lines given by equation (25) in 
the following way: If du and dv are the variations along the 
characteristic lines. 


du 

dx 



(^a* M) ) 


du 

¥ 


or (see equation (10)) 


du 

dx 


I + <*••*> I - <**>> b I 


dv 

dx 



(*-aAh ) 


av 

8 y 


(27a) 


then 


M = 8u + i Av '2 dv , a 2 8S 

dx Sx dx ay rRV an (2Tb) 


If equations (27a) and (27b) are substituted in equation (13 ), and using 
equation (26), along the characteristic line of the first family defined 

t>y 


7^ = X a = tan(p + cp) 
dx 


(28a) 


there results 


du , , dv _ a 2 dS / K 
dx ^ dx 7RV dn \H 



aw \ 1 _ 


y a 9/ 


± = 0 
H 


(28b) 


18 


MCA TN No. 1809 


and along the characteristic line of the second family defined by 



tan(cp — p) 


(28c) 


there results 


du dv a2 dS {k ^ fv clw ^ 1 

dec + dJ ~ yVR to VI ~ X *>) + + T~&i) H = 0 (28d) 


Equations (28b) and (28d) contain the term 
attack 



but at small angles of 


cos 0 = ^ cot 0 

y o9 y y 


and, therefore. 


dw 


is known at the points 


y de 

of the entropy is also known at the points 


the value of 


da 

dn 


can be determined (reference 


Pj_ and P£. The value 

and P 2 and, therefore, 
5) 


as 

dn 


AS P - AS t 


( Ip : 


- Xt 


sin 3 


cos(3 + <p) 


+ ( Ip 3 _Ip 2) 


sin 3 


cos(cp - 3) 


(29) 


P2 


From equations (28) and ( 29 ) the values of u and v can be determined 
in the first approximation for the point Pp In order to determine the 
value of w at Pp the following procedure can be used: 

If s is the projection of the streamline in the meridian plane 
considered (fig. 5) ani P 3 is a point near Py and Pp, 


dwo owe dwp 

— - = — — sin cp + ^ - cos cp 

ds by bx Y 


OWp V OWp u 

by Y bx V 
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or from equation (12) 
dw 2 

Now (fig. 5) 

% 


Ss 


yv 


?Ry 

+ 

pT 1 

11 

/dwp- 

\ . 

f y _ y- \ 

cos B 

V ds. 

' p l 

P 3 P l) 

cos(B + cp) 

+ 

11 

{ dwg s 

\ , 

: ip 3 - xp 2 ) 

cos p 

V ds. 

P 2 ( 

cos(qp — B) 


p 2 


(30) 


(31a) 


(31b) 


and 


dvjA 

dn/ E 


% " % 


/%. -y \ 

sin B 

-v \ 

sin B 


( Ip 3 p i) 

cos(<p + B) 

hd 

H 

H 

hj 

1 

hj 

cos(cp — B) 

p 2 


(32) 


Therefore, 


/8 w 2 \ 

w 2p - w 2 E + 



sin B 

( Zp 3 P 1 ) 

1 cos(qp + B) 


(33) 


The -values of u, v, w, and AS are known at the points P-^ 
and P 2 J therefore, the values of 112, v 2 , w 2 , and AS2 at the same 
points can he calculated from equations (5) and (8). Indeed, ui, V]_, 
and AS^ at those points are known from the determination of the flow 

for a = 0; therefore, from equations (30) to (33) the value of w 2 
at P3 can he determined. 

After the velocity components u, v, and w at P3 have been 
determined in the first approximation, a second approximation can he 
determined hy assuming the average values between the corresponding 
values at the points P 2 and P^ or P-j_ and P^ for all the coeffi- 
cients. After the velocity components at a point P^ have been 
obtained, the velocity components at any other point having the same x 
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and y coordinates as "but a different coordinate 9 can be 

calculated from equations ( 5 ). 

For practical calculations, equations (28) can be transformed in 
the following form; 


X a = tan(3 +- qp) 


(34a) 


^ — dcp tan 3 
V 


sin <p + 8 w 


yV be, 


tan 3 — 


sin 2 3 dS 
7 R dn 


sin 3 

cos(cp + 3 ) 


dx = 0 


(34b) 


= tan(cp - 3 ) 


(34c) 


— + dqp tan 3 
V 


tan 3 + dS 

\ y yY 80/ 7R dn 


sin 3 

cos(cp — 3 ) 


dx = 0 
(34d) 


2 _ dwg _ 1 Y2 W2 sin2p AS2 

V ds y V Yy ^ 7E y 


(34e) 


and 


where 


cp = + cq >2 cos 9 


(35) 


VO VLO 

92 = ^- COS (Pi - — Sin qp x 


(36) 


and. 


where 


V = + aV 2 cos 0 


(37) 


V 2 = U 2 cos + V 2 sin cp^_ 


(38) 


At the surface of the body the calculations are similar to the case 
of zero angle of attack because the entropy at the surface of the body 
is known in every meridian plane and the value of 9 is given. 

Equation (34e) gives the variation of W2 along the body; therefore. 
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the value of W2 can he obtained directly from another point on the 
body in the same meridian plane. 

At the surface of the shock the system of calculations is similar 
to the system for zero angle of attack. In figure 6 the point P3 is 
at the intersection of the tangents to the first characteristic surface 
at P]_ and to the shock at P2 in the meridian plane 0 = Constant. 
The equations of the shock and equation (3^b) must be verified at P^, 
which is assumed as a point of the shock in the first approximation. 

In the plane 0 = 0, w is zero and the values of V, AS, and cp 
behind the shock are functions only of the value of and for any 
value of ft, the values of Y, AS, and cp can be obtained from the 
equations of the shock 


and 


cos (n - q) = _v_ _ 
cos (ft - qp) v 0 V 0 /V 2 


(39a) 


tan(cp — a) 


?-i Ho 2 

2 M Q 2sin2(ft — a) — 1 


- 1 


tan(ft — a) 


(39*) 


AS = R ■■ log 2 ^ 7+1 

7-1 


< 7* 1 


M Q ^sin^(fi — a) 




+ 


Mo^sin^fft — a) 2 


(39c) 


IC.l + _2-Ar 


7 - 1 VL 


(39cL) 


If the plane 9 = ^ is considered, the sign of a in equations (39) 
must he reversed. 


From equations (39) the values of V and AS can he determined 
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as a function of cp; then ^ and ^ ^ as a function of cp can be 

^ dtp dqp ^ 

evaluated. Now, if cpp^ is the direction of the velocity at P-^, the 
velocity at will have the direction 


<Pp_ 


" 'Pp-L + 


Acp 


Therefore, the velocity at 
the shock of cpp^ — a and 


P^ must correspond to a deviation across 
can be expressed as 




Acp 


where Y-, is the velocity behind the shock corresponding to the 
rP-L 

direction qpp . In a similar way. 


ASp 3 




&P 


Therefore, equation (3^b) at the point P^ becomes 





(IE) - tan 




yv 


'sin cp + w ga cos 9^ sin 3 tan g 


cos(3 + cp) 


dx 







sin^ 


7 R 


!n = o 


+ 


m 
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In equation (40), Afp is the 
determined. From the value of qp 

of Yp^ can he determined; and a 


only unknown and, 
the value of ftp 


3 


therefore, can he 
and the value 


second approximation for the position 


of P-^ and its value of the velocity can he calculated if the corre- 
sponding average values between Po and P^ are assumed for cp, P, and 
all the coefficients of equation (40). 


The value of w at P can he obtained from equation(23c) in 
which T) is given by figure 6 as 


T) = sin ft 

Jp 3 


where iq = ft — a for 9 = 0 and t) = a — ft for 9 = 4.. The value 

of a corresponding to the point P^ on the shock for a = 0 is 
given by the relation 


7f V . 7p 3 


sin a sin ft 


and 


yp^ = f ^ is the curve that represents the shock for a = 0. 


PRACTICAL APPLICATION OF THE CHARACTERISTIC SYSTEM 
Graphical Numerical Method 


The analytical part of the characteristic system used for 
determining the flow field about a body of revolution at an angle of 
attack is similar to the system used for a body of revolution at zero 
angle of attack (reference 5), but the practical numerical application 
is slightly more involved. In equation (34e) the values of V 2 

and Afig must be known in order to determine the value of wq and must 
be determined from equations (37)., (36), and (8), where the values of Vq 
and z\Sq are considered known in the entire flow field and given by the 
determination from the case of zero angle of attack. In the practical 
case, however, the values of Vq, cpq, and ASj have been obtained with 
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the characteristic system only in a finite number of points at the inter- 
sections of the characteristic net; and the characteristic net for the 
case of zero angle of attack is different from the net used for the case 
of a hody with a small angle of attack. In order, therefore, to obtain 
the values of V]_ and AS]_ at the intersections of the characteristic 
lines for the case with a given angle of attack, a complicated inter- 
polation of the values Vj_ and AS]_ would he necessary if the two 

characteristic nets for zero angle of attack and for a given angle of 
attack were constructed independently. 

In order to reduce the numerical work to a minimum, the two following 
methods can he used, the first of which is practical when a graphical 
numerical calculation is performed, whereas the second can he more 
convenient when automatic computing machines are used. 

In both cases the calculations start with the determination of the 
flow at an angle of attack around a cone when the hody considered is a 
pointed-nose hody of revolution or with the determination of the shock 
at the lip of the hody if the hody is an open-nose hody of revolution. 

(The tangent to the shock at the lip can he determined with the two- 
dimensional theory.) The flow around a cone at an angle of attack has 
been determined and tabulated in reference 3 } whereas the flow for zero 
angle of attack has been tabulated in reference 6. A different method 
for determining the flow around a circular cone at an angle of attack is 
given in the appendix. It can he assumed, therefore, that the flow 
along the first characteristic line of the first family at the end of 
the conical region in the plane 9 = Constant (for example, 9 = %) 
is known (fig. 7). V 2/ 

For the practical numerical calculations a value of the angle of 
attack must he selected. In order to obtain higher precision, it is 
convenient to select a relatively high value of the angle of attack 
because in this way the differences between V and and AS and AS^ 

are large and, therefore, can he determined with sufficient precision. 

Usually, when the determination of the flow field for the case of 
zero angle of attack is made with a graphical numerical process, in 
order to avoid numerical errors of computations, the value of the 
intensity and direction of the velocity are plotted as a function of the 
position along the characteristic lines for both families of character- 
istic lines. The velocity distribution and the entropy— variation distri- 
bution along the characteristic lines and along the surface of the hody 
for the case of zero angle of attack can therefore he considered known. 

If the distribution is not given, the values of and AS^ must he 

determined as a function of x along each characteristic line of a given 
family (for example, of the second family) along the hody. 

Then the construction of the characteristic net for the selected 
angle of attack must start by drawing the first characteristic 
line P q P 2 P^ over the design of the characteristic net for zero angle 

of attack (fig„ 7). 
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From equation (34c), (34d), and (34e) the flow at Pj_ can he 

determined. From Pj_ and P2 the point P3 can he obtained in the 

first approximation as the intersection of the tangents at Pg and P^ 

to the characteristic lines. By using equations (34b), (34d), and (29), 

V, cp, and AS can he obtained in P3 as for the case of zero angle of 

attack (reference 5 ) • From the variations and al° n g "the 

line P1P3 the values of V, 9, and AS at the point P^ can be 

obtained, where P^ is obtained from the intersection of the character- 
istic line P1P3 with a characteristic line of the second family in 
the net for zero angle of attack. At the point P^, V^, AS-^, and cp-^ 

are known; and, therefore, qpg* ^2> ^2-> a* 1 *! w 2 can te obtained. From 

the values obtained from the first approximation a second approximation 
can be obtained. From Pi^ and P^ the point Pg can be determined 
in a similar way, and the flow at P7 can be calculated. Proceeding in 
a similar way, all the flow field can be analyzed. 


The equation of motion (13) can be transformed by means of 
equations (6) and (8) in a system of equations that permits a numerical 
determination of the quantities Vg > <P2-> and ^2* This system is 
numerically more involved; however, the characteristic net determined 
for zero angle of attack is used. For a small angle of attack. 


Numerical Method 


= -1_ 1 + 2 — A. (u^ + v x v 2 )a cos 9 
a 2 aj_ 2 aj_ 2 


where 



7-1 J2 

sin^p ^1 


( 41 ) 


Therefore 


— = (1 + Aa cos 9) 

a 2 a x 2 
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Substituting equations (5) in equations (7) results in the following 
expression if higher— order terms are neglected: 


Sui / _ uq 2 \ | foL f 1 V 1 2N \ u l v l + + Zl 

bx Y a l 2 / 8y Y a i 2 / ai 2 \ dy dx / y 


= — a cos 9 


^2 + ^2 
y y 


duq / 2 u 2U^ + Uq 2 A\ Su2 

dx \ aq 2 / dx 




8v n ^ 2 v 2 v i + v q^ A> 



UjVjA + UjVg + v i u 2 




u l v l / 5u 2 + 
a q 2 V 8y 8x/ 


(42) 


Because the left-hand side of equation (42) must he zero for the 
conditions at zero angle of attack. 


0U2 





uqvq /Su2 

U 



+ Q 


0 

(43) 


where 


Q = - 




-\ O 

OY-j *~ 
+ 1- 


2u 





U2V1 


u l v 2 


( 44 ) 
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or, from equations ( 10 ) and ( 36 ), 


Q = - V-I 


SYf 2ug + u^A 


6 x 


- V-, 


a l £ 


8Vi 2v 2 + VjA tp 2 V 1 8S 1 

8y a .] 2 7R dn b 


(^5) 


The value of all the coefficients at the points P-i and Pp in 

£y 

equation (43) can “be considered known “because — ^ and — ^ can “be 

ox 6 y 

considered known from the calculations for the case of zero angle of 
attack. Therefore, equation (43) can he considered an equation in which 
the characteristic lines are equal to the characteristic lines for zero 
angle of attack because the coefficients of the partial derivatives 
ClUp ClVp 


du 2 

Ik? 


dy' 8x' 


and 


dv 2 

8y 


are the same in both cases. Thus, 



X ia = tan( 9 X + p x ) 

(46a) 

X lb = tan(cp 1 - p x ) 

(46b) 


Equation (43) can be transformed by introducing the entropy 
cSSp 

gradient , and the equation of motion along each characteristic line 

on 

can be obtained. From equations ( 8 a) and (10) by means of equations (5) 
and ( 37 ), the following relation can be obtained: 


dS 

8n 


dS -i 8 Sp 

t— + — — a cos Q 

on on 


I 



SS 2 

a 

dx 


cos e 


V 1 + v 2 a 008 9 

Vi + cos 9 



+ 



\ 


cos 9 




U 1 + u 2 a cos ® 
Vj_ + V 2 a cos d 
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or, for small angles of attack. 


8S 961 8S2 

— = — — + — - a cos 9 + 

8n 8n^ 


esi 

8x 


( v 2 V 1 ~ V 1 V 2 ) 


dSn 

+ — (u2Vi - ) 


a cos 9 


J V 


(47) 


From equations (36) and (38) 


u 2 = ^2 008 9i ~ 4*2^1 sin cp^ 


(48a) 


v 2 = V2 sin cp^_ +- <p 2 v i cos 9i 


(48b) 


and 


Tjl = Y]_ sin cpj 
Ujl = cos 92 


Therefore, 


8S-, SS. 0 /8S-. u i dS. v i\ 

- ^ - nv 2 ) - ^ 0**1 - «!V 2 ) = - a * ^ Uj 


Because the term in parentheses on the right-hand side of this 
equation represents the variation of entropy along the streamline, 
which is zero, equation (47) "becomes 


as 

dn 



dS 2 

+ — ^ a cos 

onj_ 


e 


(49) 


I 
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Then, from equations (10), (41), and (49) 


(Iv ^ 8ug O' - ] SS 2 /Vg \ a^ dS^ 

8x ciy yRY^ 8n^ \Y^ J yRY^ 8n 1 


(50) 


The equations of motion along the characteristic lines defined by- 
equations (46) can be obtained by means of transformations similar to 
those of equations (27) and are 


dug + dvg + D-j. = 0 


(51) 


X la = If = tan tol + p l) (52) 

du2 + ^ la dvg + Dg dx = 0 (53) 


X lb = § = tan ^l - Pi) 


(54) 


where 



J2 A + 7 ~ 

Sni Vi_\ sii^pnLy Sn-j^ 


/ U 1 V 1 

V a i 2 - u i : 


2 + X la 


7 R 


sin^-^ 


(55a) 
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W 2 v 2 


D2 — I — + — + Q 


a i - u i 


7 - l\ 8s, 



( 5513 ) 


In equations (51) and. (53) the coefficients Dj_ and. D 2 contain the 

derivatives 1 and 1 that must he obtained from the analysis of the 
cbc dy 

case with zero angle of attack. Now, for every point P the variation 
of V-l along the characteristic line of the first family for the case 

of zero angle of attack is 


£i\ 


8v x 

dx 


+ X 


la 


la 


8 V! 

8y 


(56a) 


whereas along the characteristic line of the second family 




V^/x 


lb 


- X 8V X 

dx dy 


(56b) 


At every point P given by the intersection of two characteristic 


la 


lines X n and X-, v i n the characteristic net, the values /— i.j 
-La ID \cLx/, 

/dVn\ 

and — - 1 are known, having been obtained from the evaluation of 

V dx /X 

lb 

the following equations (reference 5): 


S = x la " tan (?l + Pi) 


(57a) 
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tan 3q 



Bin cp-) sin 3-| tan 3-] 1 
cos^ + 3 X ) y 


dSj 1 sin33i 

dn! 7R 008(92 + $2) 


(5Tb) 


= \, h = tan (92 - 32) ( 57 c) 

dx xu 

1 / dv i\ + tan 3 (' dfp l\ - I sin 9 l sin Pi tan Pi 
V2 UA lb 1 \dx/ y cos ( 9 X - 32) 

dS2 1 sin332 _ 

dn2 7R 000(92 — 32) ( 57 d) 

Therefore, the values 



can he calculated directly for every point of intersection of the charac- 
teristic line (equations (57a) and (57c)). 


After substituting the expressions of equations (48) and ( 58 ) in 
equations ( 51 ) and ( 53 ) after some simplifications and trigonometric 
transformations, the following equations can be obtained: 


A, 


la 


= ^ = tan (<Pi + p l) 


(59a) 
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dY2 

Y^- ~ d( P 2 tan 


+ 


AlPl 




sin 


9 X + 


cos 



dx 


+ 


T dx 

sin 3-^ cos 2 3 1 cos^ 3 1 + cp-jJ 


= 0 


(59b) 


x lb ~ = tan ( < Pi “ 


(59c) 


dYo ApPp Ao /V Q w n \ 

+ dqpp tan 3 n — dx ^ f sin cp, + © 0 cos qp, + - 2 . | dx 

Yi d 1 y \V X 1 Y2 Y 1 ) 


T dx 


sin 3 1 003^3^^ cos^cp^ — 3^) 


= 0 


(59d) 


where 


A _ tan 3 ^ sin 3 l 
1 cos^-l + 3 X ) 


( 60 a) 


_ tan 3i sin 3 ] 
cos^ - 3 x j 


( 60 h) 


Ss P 

P-L = — — cos 3 t sin 3 , — 
Sn-i x 


^1. 


i 1 * ^ l )°° 3 ^ 8ln P 1 - ^ 


dSj 

dn n 


( 60 c) 


SSp 

P 0 = — - cos 3 t sin 3 , - 

on-j^ 


Vi + 


7 - 1' 

sin : 


- — ] cos 3 , sin 3 , + qpp 

2 Pl) 1 1 2 


dS-, 


<3n n 


( 60 d) 
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cos^ + <p lrf ) 


Wl + . 7-1 

v x l 2 Sln2p-| 


jsin P-j. + cp 2 cos P-l 



008(^2 — 



a^inL )° ln " ”2 008 p i 


( 60 e) 


dS-, 

The coefficients A-, and Ao and the value of — =• have been 

1 d dn-L 

determined for the flow at zero angle of attack and 



The practical use of equations (59) is identical to the use of the 
corresponding equations (57) for the case of zero angle of attack. 

(See reference k.) 


CONCLUDING REMARKS 


The method of characteristics has been applied to bodies of 
revolution at a small angle of attack. Only the first-order effects of 
the angle of attack have been considered. The system developed takes 
into account the effects of the entropy variations on the flow phenomena 
and determines a flow that exactly satisfies the boundary conditions 
within the limits of the simplifications assumed. 

The application of the method to practical problems has been 
discussed and two systems are given. The first method is numerical and 
analytical and requires less numerical computation but requires the 
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construction of another characteristic net, whereas the second method is 
only numerical and uses the characteristic net and some of the numerical 
computations made for the calculations for zero angle of attack. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va„, November 22, 19^8 
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APPENDIX 

DETERMINATION OF FLOW PROPERTIES AROUND A CIRCULAR CONE 
AT A SMALL ANGLE OF ATTACK 


Assume a polar coordinate system r,\|r,0. Call v r the velocity 
in radial direction, v n the velocity in normal direction to r in 
the meridian plane 9 = Constant, and w the component normal to the 
meridian plane (fig. 8); that is. 


v r = xr 


dr 

dt 


v n = 


v = 


r d^ 
dt 

r d9 sin t|t 
dt 


If the phenomenon is conical. 


and 


or 

= 0 

Sr 

^-0 


8p 

cir 


= 0 


= 0 

8r 


Therefore, Euler’s equations are 


v n 8v r 


w 


OVy 


'n 




6 * 


r sin 


89 


= 0 


(62a) 


la ^ 


W 


S^ r sin o9 pr S^ 


8la + _L & + irZn ~ y2cot ^ = 0 


(62b) 
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V n dw W dw 1 dp V.„W + V-W cot \|r 

+ — + — + -£ s = o 


r dv r sin ijr bd rp sin \|/ d9 


(6 2c) 


The continuity equation is 


2pv sin i|r + v sin Tjf ^- + p sin t + v n p cos + +■ ¥ ^ + ^ 
n St d* n d<9 d 9 


= 0 


(63) 


and the energy equation is 


-.7 A °P _I^£! 
7 — 1 \p d9 p^ d0 


m -L '^X. + „ 5Xn + w 3* 

V r 89 n 89 89 


(64a) 


7 /l cip P d P \ _ ( 6v r 6v n dw 


= — v + T „ — + w 


7 — 1 l p di|r p2 di|r/ \ r di|f n dty di|f 


(64b) 


Combining equations (62) to (64) results in 


T - ia 


a^ / sin <30 ' 


<3w / ___ 


2w n dv n + d^ _ ^ _ w cot +] = 0 


a2 sin If d 9 a2 Is in 4' bd dty 


(65) 


For small angles of attack the velocity components can he expressed 
in the form (references 1 and 3) 


v = v + av cos 0 
r r x r 2 


v n = % + av n 2 cos 9 


(66) 


w = aw2 sin 6 
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when the second and higher order terms of the angle of attack have been 
neglected. Equation (65) at small angles of attack becomes 




+ v. 


n 


8v 

cot ^ + — — 
8^ 


v n 2 j , dw 
a^ / sin tJt 89 


0 


(67) 


At small angles of attack. 


p = p-^ + ap 2 cos d 


p = + a,p 2 cos 9 


The shock is a circular conical shock having its axis inclined at an 
angle t) with the axis of the cone. The quantities with subscript 1 
are the quantities corresponding to the case of zero angle of attack, 
(indeed, the cone is a particular body of revolution; and, therefore, 
the considerations made for the case of bodies of revolution are still 
valid. ) 

From equations (62a), (62b), and (64b), there results at small 
angles of attack 


1 d£ _ z = 0 

p p 8 i 

In the meridian plane 9 = Constant, therefore, the transformation 
behind the shock is isentropic for small angles of attack. If a AS 2 cos 9 
is the variation of entropy in a direction normal to the meridian 
plane 9 = Constant, 


8S 

r sin \|r 8© 


a AS2 cos ^ 


Because the shock at small angles of attack has circular cross section. 


7-1 


as 2 



E 


= Constant 
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and from equations (62c) and (64a) 



dwo 

V n 3in * + w 2 


( v r sin * + v n cos ^ + v r v r 2 + 


( 68 ) 

If the variation of entropy is small and the term AS 2 can be neglected, 
equation (68) becomes 


Wo sin Mr = — v 
2 r 2 


( 69 ) 


Equation (6 7) can be written in the following form: 



+ 



— v r - v n cot i|r 


aw 2 cos 9 
sin i|r 


(70) 


By use of equation ( 67 ) and by considering the conditions for zero angle 
of attack, equation (70) becomes: 



= —v 


"2 


cot i|r 




Y r ± + v n x cot ♦ 





( 71 ) 


NACA TN No. 1809 


39 


Equations (68) , (70), and (71) permit the determination of the flow 
around the cone at an angle of attack by means of a step— by— step calcu- 
lation when the calculation for a = 0 has been performed. Consider 
the hodograph plane uv, and consider the variation of velocity 
components v r and v n in a meridian plane 9 = Constant (fig. 9)* 

Assume that for a given value of t|r a and 9 & the velocity 
components v r , v n , and w are known. Point P a of the hodograph 
diagram represents the velocity vector 0P a corresponding to the velocity 

at every point of the space of coordinate t a in the plane 6 a = Constant 
whereas 0Q a represents the values of ( v n)p and - Qa^a represents the 

values of (v r ^p . Now , the radius of curvature R a of the hodograph 
diagram is along the line Q a P a and has a value given by (reference 1) 



and, therefore, can be determined from equation (70). At P Q the values 
of v and v_ are known from the calculation for a = 0; 
r l n l 

therefore, av and av can be determined from equations (66). 
r 2 **2 

Equation (71) can be used in place of equation (70) in the following way: 
The vectors 0Q a and Q a P & represent the values of v^ and v r ^ 

at P : the vector 0P Q in the hodograph diagram gives the values 

a ’ a. 

of Vg j and, therefore. 



can be obtained from equation (71)* 

Now at any point P a the radius R a , given from equation (70), or 

the radius R 2 > given from equation (71), is known; therefore, from the 
a 

quantities at P a the quantities at Pb of coordinate ^ = t a + At 
can be obtained by constructing a circle of center C a (where C & P a = R a 
or P 2 ^ through the point P a until the point P^, along the line C, 
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which is a straight line from C & , is inclined by i|r + A\|r with the 
u-axis. Therefore,, 


C v n)^ +A \|r = C v n)^ cos At) + (R - v r )^ sin At) 
( v r)*+ Ailr = ( v n)y 8ln AT) “ (R ” v r)i|r coS At) + 


Inasmuch as the values of v r and v n at P b have been obtained, 

the values of v and v can be determined by differences from the 
r 2 "2 

values for a = 0 with the use of equations (66). (if equation (71) is 
used, the values of v and v n are obtained directly. \ With 
, .. _ _ 6w2 / 


equation (68) the value of 


can be calculated at P D , and the value 


of W2 at P^ can be obtained. Indeed, AS2 is constant and has been 


determined from the conditions at the shock. In a similar way, all the 
holograph diagram can be constructed. If necessary, for every point P-^ 

a second approximation can be determined. 


The calculation of all the flow field must start at the shock. For 
the calculations it is convenient to choose a coordinate system having 
the axis of the conical shock as the axis of polar coordinates. In this 
case, the velocity components v r , v n , and w behind the shock can still 

be expressed in the form of equation (66). Indeed, from equations (15) 
and (19), 



v_ + a cos 
n l 




y + ay cos 
n lg n2 


\ 


(72a) 



a cos 




v 2 r 


+ ay 


cos 9 


(72b) 
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w = a sin 9 ! w 2 


L n - -±— a 


V 2 sin * a n l 


= aw 9 sin 0 


(72c) 


where v„ , v_ , and. w„ are the components referred to the axis of the 
n s r s ° 

conical shock, whereas the components v.^, v^, and w^ are referred 
to the axis of the body. Indeed, r)/a is constant. 

The calculations start at the shock. After determining the flow 
field for zero angle of attack, the angle of the conical shock is 

known and the velocity components v n and v r with respect to the axis 

of the shock for every value of i are also known. In order to determine 
the flow for the case of a small angle of attack, the direction of the 
undisturbed velocity must he rotated at a small angle a — T| with respect 
to the axis of the shock (fig. 10). The value assumed for a - ti fixes 
the value of a for which the calculations are performed. (This value 
of a is not yet known hut is obtained as a result of the calculation.) 

For the value of a — t] chosen, the components v r and w s behind 

s 

the shock can he determined from equation (l4b) ^v r = v T ’j and 
equation (14c) (w» = w s ) ; whereas v n ^ can he determined from equa- 

tion (21a) ^v^ = - vn") and V N* is g iven ^ equation (14a). The 

value of entropy a AS 2 cos 9 can also he determined from the equation 
of the shock, for example, from the difference between AS and ASp. 

When v r , v n , w 2 , and AS 2 are known behind the shock, all the flow 

field can he obtained by means of equations (68) and (70) or (71). The 
holograph diagram can he constructed, for example, in the plane 9 = 


The axis u has been chosen in the direction of the undisturbed 
velocity for zero angle of attack that corresponds to the axis of the 
shock for a = 0. For a the undisturbed velocity has been rotated 
at a - t) with respect to the u-axis (fig. 10); therefore, the axis of 
the shock has not been changed. The velocity 0P a behind the shock of 
figure 9 must be decomposed (l) in a component P a Q a inclined at V 3 

corresponding to v r if equation (70) i s used or to v r ^ if 
equation (71) is used and (2) in a component 0 a Q a corresponding to v^ 


k2 


MCA TN No. 1809 


In this way, the values of 


'n 


and v_ 
n i 


that must he used, in 


equations (68), (70), and (71) are the values obtained from the calcu- 
lations for zero angle of attack at the same value of \|r 
(jthat is, ^ = i|r s for P a ) . Because the calculations start at the 
shock, the construction of the holograph diagram must he performed in 
the direction of decreasing values of \(r . At the surface of the body 
for 9 = A., the component v n must he zero; therefore, when the radius 

of the holograph diagram passes at the origin of coordinates u and v, 
the corresponding value of i|r is equal to i|r 0 + q where \[r 0 is the 
angle of the cons (fig. 10). Because \|r 0 is known, the value of r| 
and, therefore, of a can he determined. 


The components v r and v n in the plane 0=0 or 


0 = - do not 
2 


change when, for the axis of reference, the axis of the body is assumed; 
hut the corresponding value of t|t is increased at q (fig. 10). The 
value of Wg changes; the value of w 2-^" can determined from the 

value of wg by means of equation ( 72 c). 


For practical calculations it is convenient to use nondimensional 
coefficients obtained by dividing all the velocity components by the 
limiting velocity . The expression a/V ^ can be obtained from 
equation (39d). 


For small values of a, the values of 


Vr 2* Vn 2* w 2> 


and 


q/a are 


independent of a; and, therefore, the flow for every other value of a 
can be obtained from this determination. The calculations can be 
graphical or analytical. 
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Figure 1.— Cylindrical coordinates and corresponding velocity components. 
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Figure 2.— The shock surface at a = 0 


and a ^ 0. 
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Figure 3-— The velocity components in front of and behind the shock. 
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Figure 4.- The velocity components behind the shock for a ^ 0 as a 
function of the velocity components behind the shock for a = 0. 
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Figure 5 .- The analysis of the flow with the characteristics system. 
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Figure 6.— Determination of a point of the shock. 
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Characteristic net for 
a small angle of attack 

Characteristic net 
for zero angle of attack 


e = tt 



Figure 7-— Scheme of the characteristic net. 
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Figure 8.— Polar coordinates and velocity components. 
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Figure 9-— The hodograph diagram. 
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Figure 10.- Choice of the axis for the holograph diagram. 


